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Abstract 

Many-body quantum-mechanical scattering problem is solved asymptotically when 
the size of the scatterers (inhomogeneities) tends to zero and their number tends to 
infinity. 

A method is given for calculation of the number of small inhomogeneities per 
unit volume and their intensities such that embedding of these inhomogeneities in a 
bounded region results in creating a new system, described by a desired potential. 
The governing equation for this system is a non-relativistic Schrodinger's equation 
described by a desired potential. 

Similar ideas were developed by the author for acoustic and electromagnetic (EM) 
wave scattering problems. 

PACS 43.20.+g, 62.40.+d, 78.20.-e. 
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1 Introduction 

In a series of papers, cited in the references, and in the monograph [8] the author has 
developed wave scattering theory by many small bodies of arbitrary shapes and, on this 
basis, proposed a method for creating materials with some desired properties, in particular, 
with a desired refraction coefficient. The goal of this paper is to describe a new method for 
creating materials with a desired refraction coefficient and to compare it with the method 
proposed earlier. The new method is applicable to the problems which deal with vector 
fields and tensorial refraction coefficients. Our presentation deals with the simpler case of 
scalar fields and uses paper [27] . 

Let us formulate the statement of the problem. Let D C i 3 be a bounded domain 
filled in by a material with a known refraction coefficient. The scalar wave scattering 
problem consists of finding the solution to the Helmholtz equation: 

L u := [V 2 + k 2 nl(x)}u = in R 3 , Lung (a) > 0, (1) 
u = u + v, u :=e ika - x , (2) 
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e 1 x 

v = Aq{P,ol) h o(-), r := \x\ — > oo, (3 := —. 



(3) 




The function tiq(x) is assumed bounded, Riemann-integrable, 



ng(s) = l in D':=K 3 \D. Imra§(x) > 0. 



(4) 



The function ^4(/3,a) is called the scattering amplitude. The wavenumber fc = where 
A is the wavelength in .D', a £ S* 2 is the direction of the incident plane wave uq, S 2 is 
the unit sphere in IR 3 , (3 6 5* 2 is the direction of the scattered wave. The solution to 
problem {I])-© is called the scattering solution. It is well-known ([5]) that this solution 
exists and is unique under our assumptions. These assumptions are: n^{x) is bounded: 
su Pa;eiR 3 l n o( x )l — n o = const, Imno(x) > 0, n^{x) = 1 in D' . Since the solution u is 
unique, the corresponding scattering amplitude Ao((3,a) is determined uniquely if ng(x) 
is given. We assume k > fixed and do not show the dependence of Aq(/3,o) on k. 
The operator Lq at a fixed k > can be considered as a Schrodinger operator Lq = 
V 2 + k 2 - qo(x), where q (x) := k 2 - k 2 nl(x). 
Problem RC: 

We want to construct a material in D with a desired refraction coefficient, that is, with 
a desired function n 2 (x), Imn 2 (a;) > 0. 

Here RC stands for refraction coefficient. 

Why is this problem of practical interest? 
We give two reasons. 

First, creating a refraction coefficient such that the corresponding material has negative 
refraction is of practical interest. 

One says that a material has negative refraction if the group velocity in this material 
is directed opposite to the phase velocity. 

Secondly, It is of practical interest to create a refraction coefficient such that the 
corresponding scattering amplitude A(f3) = A{j3, a) at a fixed a G S 2 (and a fixed k > 0) 
approximates an arbitrary given function /(/?) 6 L 2 (S 2 ) with any desired accuracy e > 0. 

This problem we call the problem of creating material with a desired wave-focusing 
property. 

In Section 2 we explain how to create a material with a desired refraction coefficient. In 
Section 3 we compare our recipe for creating material with a desired refraction coefficient 
with the recipe given in [24J. 

2 Creating material with a desired refraction coefficient by 
embedding small inhomogeneities into a given material 

The basic idea of our method is to embed in D many small inhomogeneities q m (x) in such 
a way that the resulting medium can be described by a desired potential q(x). 
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Let us assume that Pm(%) is a real- valued compactly supported bounded function, 
which is a sum of small inhomogeneities: 



M 

Pm = ^2 9m(x), 

m=l 

where q m (x) vanishes outside the ball B m := {x : \x — x m \ < a} and q m = A m inside B m , 
1 < m < M, M = M(o). 
The question is: 

Under what conditions the field um, which solves the Schrodinger equation with the 
potential pm{x), has a limit u e (x) as a — > 0, and this limit u e (x) solves the Schrodinger 
equation with a desired potential q(x) ? 

We give a complete answer to this question in Theorem 1 below. 

The class of potentials q, that can be obtained by our method, consists of bounded, 
compactly supported, Riemann-integrable functions. It is known that the set of Riemann- 
integrable functions is precisely the set of almost everywhere continuous functions, that 
is, the set of bounded functions with the set of discontinuities of Lebesgue measure zero 
in M 3 . These assumptions on q are not repeated but are always valid when we write "an 
arbitrary potential" . 

In fact, a more general set of potentials can be constructed by our method. It is 
mentioned below that for some class of unbounded potentials, having local singularities, 
which are absolutely integrable, our theory remains valid. 

Our result is as follows: 

Assume that q(x) is an arbitrary Riemann-integrable in D potential, vanishing outside 
D, where D is an arbitrary large but finite domain, and the functions A{x) and N(x), 
which we can choose as we wish, are chosen so that A(x m ) = A m and A(x)N(x) = q(x), 
where N(x) > 0. 

Then the limit u e (x) of um(x) as a — >• does exist, and solves the scattering problem 
with the desired refraction coefficient n 2 (x): 

V 2 u e + k 2 n 2 (x)u e := V 2 u e + k 2 u e — q(x)u e = 0, (5) 

u = u + v, (6) 

where uq solves problem ([5]) and v satisfies the radiation condition. 

The notation u e (x) stands for the effective field, which is the limiting field in the 
medium as M — > oo, or, equivalently, a — > 0. Under our assumptions (see Lemma 1 
below) one has M = 0(\). 

The field um is the unique solution to the integral equation: 

r e ik\x-y\ 

u M {x) = u (x) - \ / g{x,y,k)q m {y)u M {y)dy, g(x,y,k) = — -, (7) 

' d 4tt\x — y\ 



m=l 



where uq(x) is the incident field, which one may take as the plane wave, for example, 
uq = e tka ' x , where a £ S 2 is the direction of the propagation of the incident wave. 
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We assume that the scatterers are small in the sense ka « 1. Parameter k > is 
assumed fixed, so the limits below are designated as limits a — > 0, and condition ka « 1 
is valid as a — > 0. 

If ka « 1, then the following transformation of ([7]) is valid: 

u M (x) = uq(x) - 22 — J A m u M {x m ) / I _ I [1 + o(l)\. (8) 

m= l 47F J\y-x m \<a \ x V\ 

To get ([SD we have used the folowing estimates: 

\x — x m \ — a <\x — y\ <\x — x m \ + a, \y — x m \ < a. 
These estimates imply 

e ik\x- y \ = e «l*-*ml[i + (i)] 9 

provided that \y — x m \ < a and a — > 0. 

We also have taken into account that, as a — > 0, one has: 

max \u M (x) - u M {xm) \ = o(l). 

A proof of this statement is given in [12j. 

Here is a sketch of another proof of the above statement. Equation (|7|) has a unique 
solution because it is an equation with a compactly supported bounded uniformly with 
respect to M potential p. Its solution satisfies a homogeneous Schrodinger equation and 
the radiation condition. Thus, this solution is unique. This solution is uniformly (with 
respect to M) continuous in D. 

Yet another argument, different from the ones, given above, can be outlined as follows. 
The limiting funtion u e (x) is in H 2 (D) by the standard elliptic regularity results, so it is 
continuous in D. The function um converges to u e uniformly. Therefore, it satisfies the 
above inequality as M — >• oo, or, equivalently, as a — > 0. 

The fourth proof of the same statement can be based on a result from |26j . Namely, 
a justification of a collocation method for solving equation (17) (see below) for the lim- 
iting field u e is given in [26j. From the arguments, given there, one obtains the uniform 
convergence of um to u e : 

lim max \um(x) — u e (x)\ = 0. 
Let us return to equation ([8|). 

We want to prove that the sum in ([8]) has a limit as a —> 0, and to calculate this limit 
assuming that the distribution of small inhomogeneities or, equivalently, the points x m , is 
given by the following formula: 

M(A) = \V(a)\- 1 [ N(x)dx[l + o(l)] a^0, (9) 

J A 

where V(a) = 4vra 3 /3, N(x) > is an arbitrary given function, continuous in the closure 
of D, < N(x) < p < 1, p is a constant arising in the problem of spheres packing, and 
A/"(A) is the number of small inhomogeneities in an arbitrary open subset A C D. 
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The maximal value of the constant constant p is the maximal ratio of the total volume 
of the packed spheres divided by |A|, where |A| is the volume of A. The total volume of 
the balls, embedded in the domain A, is equal to 

V A := V(a)Af(A) = [ N{x)dx[l + o(l)]. 

J A 

If diam(A) is sufficiently small and x € A is some point in A, then 

V A = N(x)\A\[l + o(l)] < p|A|. 
Since the domain A is arbitrary, one lets a —¥ and concludes that 

N(x) < p < 1. 

It is conjectured that maximal value of p is equal approximately 0.74. (see, e.g., |32j). 

There is a large literature on optimal packing of spheres (see, e.g., [3], [32]). For us 
the maximal value of p is not important. What is important is the following conclusion: 

One can choose N(x) > as small as one wishes, and still create any desired potential 
q(x) by choosing suitable A(x) > 0. 

The total number M = M{D) of the embedded small inhomogeneities is M = 0(\), 
according to ([9]). 

Our basic new tool is the following lemma. 

Lemma 1. If the points distributed in a bounded domain D CM 3 according to 

([9]), and f{x) is an arbitrary Riemann-integrable in D function, then the following limit 
exists: 

M 

lim V f(x m )V(a) = / f(x)N{x)dx. (10) 

The class of Riemann-integrable functions is precisely the class of bounded almost 
everywhere continuous functions, i.e., bounded functions with the set of discontinuities of 
Lebesgue measure zero in M 3 . The class of such functions is precisely the class of functions 
for which the Riemannian sums converge to the integrals of the functions. 

The result in (jlOp we generalize to the class of functions for which the integral in the 
right-hand side of (|10p exists as an improper integral. In this case f(x) may be unbounded 
at some points y, but the limit lim^o j£> f(x)N(x)dx exists and 

lim / f{x)N{x)dx := / f(x)N(x)dx, 

where D$ := D\ B(y, 5), and B(y, 5) is the ball centered at y E D and of radius 5. In this 
case the sum in (jlOh is defined as follows: 

M 

lim V f(x m )V(a) := lim lim V f{x m )V(a). (11) 

a->0 L — ' 5^0 a^O ^-^ 

m=l x m £D s 
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The same definition is valid for the conclusion of Theorem 1, which is our result. 

Theorem 1. If the small inhomogeneities are distributed so that ([9]) holds, and 
q m (x) = if x B m , q m {x) = A m if x £ B m , where B m = {x : \x - x m \ < a}, 
A m := A(x m ), and A(x) is a given in D function such that the function q{x) := A(x)N{x) 
is Riemann-integrable, then the limit 

lim um (x) = u e (x) (12) 

a— s>0 

does exist and solves problem dS])-© with the potential 

q(x) = A(x)N(x). (13) 



There is a large literature on wave scattering by small inhomogeneities. A recent paper, 
which we use, is [27]. Some of the ideas of this approach were earlier applied by the author 
to scattering by small particles embedded in an inhomogeneous medium (see [3]- [23]). 

We apply Lemma 1 to the sum in ([8]), in which we choose A m := A(x m ), where A(x) 
is an arbitrary continuous in D function which we may choose as we wish. A simple 
calculation yields the following formulas: 

/ \x - y^dy = V(a)\x - x m \~ l , \x - x m \ > a, (14) 

J \y-x m \<a 



and 

\x — y\~ 1 dy = 27r(a 2 — — ), \x — x m \ < a. (15) 



\x - |2 



\y-x m \<a 3 



Therefore, the sum in ([8]) is of the form (|10p with 

f( X m) = -T-, rA(x m )u M (x m )[l + o(l)]. (16) 

4:7T\X — X m \ 

Applying Lemma 1, one concludes that the limit u e (x) in (|12p does exist and solves the 
integral equation 

e ik\x-y\ 



u e (x)=u (x)- — rq(y)u e (y)dy, (17) 

Jd 4vr|x - y\ 

where q(x) is defined by formula (|13j) . 

Applying the operator V 2 + k 2 to dUJ, one verifies that the function u e {x) solves 
problem ©-(H]). 

Thus, the conclusion of Theorem 1 is established. 

In (fl~6|) f(x m ) depends on a through UM{x m )- Nevertheless Lemma 1 is applicable 
because um{x) converges to the solution u e (x) of eqution (fT7|) as a — > 0. This follows from 
the results in |26] and [12]. These results give two independent ways to prove convergence 
of um to u e as a — > 0. The results in [12] yield compactness of the set {um{x)}, while 
the results in |26] convergence of a version of the collocation method for solving equation 
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(|17|) is proved. This version of the collocation method leads to solving the following linear 
algebraic system 

n 

uj = u j - ^2 9jmqmU m V(a), 1 < j < m, (18) 

ik\xj— x m \ 

where Uj := u{Xj), q m := q(x m ), and gj m := %^ x ._ x y It is proved in [26] that the linear 
algebraic system (|18|) has a solution and this solution is unique for all sufficiently small a, 
and the function 

n 

u {n \x) := UjXj(x) 

m=l 

converges to u e (x) uniformly: 

lim \\u e {x) -n (n) (x)|| c(D) = 0, 

where Xj( x ) = 1 x G D, Xj(x) = x £ D. 

In particular, it is proved in [26j that a collocation method for solving equation (|17p 
leads to sums, similar to the sum on the left-hand side of (|10p . and the collocation method, 
studied in [26], converges to the unique solution of the equation (fTTj) . This means that 
the function um(x) converges to u e (x). 

Remark 1. The quantities A m in the definition of q m can be it tensors. In this case 
the refraction coefficient n 2 (x) of the created material is a tensor, and u is a vector field. 

Remark 2. The recipe for creating materials with the desired refraction coefficient by 
the method, based on Theorem 1, can be summarized as follows: 

Recipe 1: Given a bounded domain filled with the material with the refraction coefficient 
71q(x) ; one embeds the small inhomogeneities q m , q m = for \x — x m \ > a, q m = A m := 
A(x m ), where the points distributed by the rule ([9]), and one chooses A(x) and 

N(x) > 0, so that A(x)N(x) = q(x), where q(x) is the desired potential, or which is the 
same, the desired refraction coefficient n 2 (x) := 1 — k~ 2 q(x). The material, created by the 
embedding of these small inhomogeneities, in the limit a —¥ 0, has the refraction coefficient 
n 2 (x) = 1 - k~ 2 q(x). 

One can fix N(x) > and choose A(x) such that the function q(x) = A(x)N(x) is a 
desired function, so that n 2 (x) = 1 — k~ 2 q{x) is a desired function. 
We call this recipe Recipe i, where i stands for inhomogeneities. 

3 A discussion of two recipes 

In [24] a different recipe was proposed for creating materials with a desired refraction 
coefficient. We call the recipe from [24] Recipe p, where p stands for particles. Recipe 
p consists of embedding small particles, for example, balls of radius a, ka « 1, into a 
given material according to the distribution law, similar to the one in ([9]), but with a 2 ~ K 
in place of V(a), where k € (0,1] is a parameter experimenter can choose at will. The 
physical properties of the embedded particles are described by the boundary impedance 
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Cm = h a™ ■ The function h(x), determining the boundary impedances, can also be chosen 
as one wishes. We refer the reader to the paper |24] for details. Numerical results on the 
implementation of Recipe p are given in [2] and [1]. 

There are two technological problems that should be solved in order that Recipe p be 
implemented practically. 

Technological Problem 1: How does one embed many small particles in a given 
material so that the desired distribution law is satisfied? 

Technological Problem 2: a) How does one prepare a small particle with the desired 
boundary impedance ( m = h ^™^ ? 

b) How can one prepare a small particle with the desired dispersion of the boundary 
impedance, that is, the desired w- dependence of h = h(x,uj)? 

The motivation for the last question is the following: by creating material with the 
refraction coefficient depending on the frequency u in a desired way, one can, for instance, 
prepare materials with negative refraction. 

The first problem, possibly, can be solved by stereolitography. 

The second problem one should be able to solve because the limiting cases Cm = 
(hard particles) and Cm = oo (soft particles) can be prepared, so that any intermediate 
value of Cm one should be able to prepare as well. 

The author formulates the above technological problems in the hope that engineers 
get interested and solve them practically. 

In Recipe p the number M of the embedded particles is M := Mi = 0( ^i K ), while in 
Recipe i the number M is M := M\ = O(-^j). Thus, M% » Mi as a — > 0. So, in recipe 
2 one has to embed much smaller number of small particles than in Recipe i. This is an 
advantage of Recipe p over Recipe i. The disadvantage of Recipe p, compared with Recipe 
i, is in the possible practical difficulties of creating particles with the prescribed very large 
boundary impedances Cm = O(A-) when a is very small. 

3.1 Negative refraction 

Material with negative refraction is, by definition, a material in which group velocity is 
directed opposite to the phase velocity, (see [1] and references therein). Group velocity is 
defined by the formula v g = VkC<j(k). Phase velocity v p is directed along the wave vector 
k° = to. In an isotropic material u = w(|k|), and um(x, u) = c|k|. Differentiating this 
equation yields 

V k a;[n(x,u;) + = ck°. (19) 

Thus, 

v a [n + (20) 

Wave speed in the material is |v p | = , where c is the wave speed in vacuum (in D'), 

and n(x,uj) is a scalar, and v p is directed along k°. 

For v g to be directed opposite to v p it is necessary and sufficient that 

dn 

n + u—<0. (21) 

OOJ 
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If the new material has the function n{x,oj) satisfying (|2ip . then the new material has 
negative refraction. 

One can create material with the refraction coefficient n 2 (x,uj) having a desired dis- 
persion, i.e., a desired frequency dependence, by choosing function h = h(x,u)) properly. 

Let us formulate the technological problem solving of which allows one to implement 
practically Recipe p for creating materials with negative refraction. 

Technological Problem 3: How does one prepare a small particle with the desired 
frequency dependence of the boundary impedance Cm = -^-^r-^-, where h(x,uj) is a given 
function? 

This problem is the same as Technological Problem 2b), formulated above. 
3.2 Wave- focusing property 

Let us formulate the problem of preparing a material with a desired wave-focusing property. 
The refraction coefficient is related to the potential by the formula 

q(x) = k 2 — k 2 n 2 (x). 

We assume in this Section that k > is fixed, and a € S 2 , the incident direction (that is, 
the direction a of the incident plane wave e tka ' x ) is also fixed. 

Problem: Given an arbitrary fixed f(/3) € L 2 (S 2 ), an arbitrary small fixed e > 0, an 
arbitrary fixed k > 0, and an arbitrary fixed a S S 2 , can one find q € L 2 (D), q = in 
D' = R 3 \D, such that 

K(/3)-/(/3)|| L2(s2) <e, (22) 

where A q (/3) := A q (J3,a,k) is the scattering amplitude, corresponding to q, at fixed a and 
k, and q := k 2 — k 2 n 2 (x) ? 

The answer is yes, and an algorithm for finding such a q, and, therefore, such an n 2 (x), 
is given in |17j . See also [18], |20| . The above problem is the inverse scattering problem 
with the scattering data given at a fixed k and a fixed direction a of the incident plane 
wave. 

There are many potentials q which solve the above problem. It is possible to choose 
from the set of these potentials the one which satisfies some additional properties, for 
example, one can choose q to be arbitrarily smooth in D, one can try to choose q with 
non-negative imaginary part, etc. 

The desired radiation pattern /(/?) can, for example, be equal to 1 in a given solid 
angle and equal to outside this solid angle. In this case the scattered field is scattered 
predominantly into the desired solid angle. 

This is why we call such a material a material with wave-focusing property. 
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